High incidence of imperforate vagina in ADGRA3-deficient mice

Background Ten percent of the female population suffers from congenital abnormalities of the vagina, uterus, or oviducts, with severe consequences for reproductive and psychological health. Yet, the underlying causes of most of these malformations remain largely unknown. ADGRA3 (GPR125) is involved in WNT signaling and planar cell polarity, mechanisms vital to female reproductive tract development. Although ADGRA3 is a well-established spermatogonial stem cell marker, its role within the female urogenital system remains unclear. Results In this study, we found Adgra3 to be expressed throughout the murine female urogenital system, with higher expression pre-puberty than after sexual maturation. We generated a global Adgra3−/− mouse line and observed imperforate vagina in 44% of Adgra3−/− females, resulting in distension of the reproductive tract and infertility. Ovarian morphology, plasma estradiol, ovarian Cyp19a1, and vaginal estrogen receptor α (Esr1) expression were unaffected. However, compared to controls, a significantly lower bone mineral density was found in Adgra3−/− mice. Whereas vaginal opening in mice is an estrogen-dependent process, 17β-estradiol treatment failed to induce vaginal canalization in Adgra3−/− mice. Furthermore, a marked reduction in vaginal and ovarian progesterone receptor expression was observed concomitant with an upregulation of apoptotic regulators Bcl2, Bid, and Bmf in adult Adgra3−/− females with a closed vagina. Conclusions Our collective results shed new insights into the complex mechanisms by which the adhesion receptor ADGRA3 regulates distal vaginal tissue remodeling during vaginal canalization via altered sex hormone responsiveness and balance in apoptotic regulators. This highlights the potential of ADGRA3 as a target in diagnostic screening and/or therapy for obstructive vaginal malformations in humans. Supplementary Information The online version contains supplementary material available at 10.1186/s12915-024-01873-6.


Background
Normal development and function of the female reproductive tract are indispensable for mating, fertilization, implantation, embryonic development, and safe delivery of offspring.However, up to 10% of women suffer from congenital abnormalities of the vagina, uterus, or oviducts [1][2][3], which has devastating consequences for their reproductive and psychological health [4].The most common congenital obstructive anomalies of the reproductive tract are imperforate hymen (observed in 1:1000 females), followed by a transverse vaginal septum and vaginal atresia [5,6].Despite occurring at high frequencies, the etiologies of most vaginal abnormalities are unknown [7].An increased effort toward understanding basic biology and underlying causes of vaginal malformations is necessary to improve diagnosis and appropriate treatment.
Although the anatomy of mice differs slightly from that of humans, the signaling pathways are largely conserved, making mouse models an invaluable tool for investigating female reproductive tract development [8,9].The female reproductive tract originates from the mesodermal Müllerian ducts and is guided by the Wolffian ducts to grow caudally, eventually reaching the endodermal urogenital sinus [9].Fusion of the Müllerian ducts results in the formation of the uterovaginal canal.The urogenital sinus subsequently forms the sinovaginal bulbs that develop into the solid vaginal plate after proliferation.Growth and canalization of the vaginal plate occur simultaneously to form the lower part of the proximal vagina and vaginal orifice [9].In mice, this process occurs at puberty via estrogen-mediated induction of apoptosis of the distal vaginal epithelium to generate the vaginal entrance [10][11][12].Thus, the opening of the murine vagina remains closed at birth and is sensitive to local and systemic alterations in estrogen signaling that accompany sexual maturation [12][13][14].The exact apoptotic mechanism that produces the vaginal opening remains unclear, but Müllerian duct development depends on proper WNT signaling [9,15].In support of a complex mechanism of vaginal development, obstructive vaginal malformations and failure of opening have been reported in several mouse strains, including those with genetic mutations in components of the WNT pathway cascades [16][17][18][19][20][21].They have also been reported in mice with mutations in tumor suppressors (LHFPL2 G102E and trBraca1 + p53 +/− ) [22,23], overexpression of apoptotic regulator BCL2 [11], and mono, dual, or triple knockout of other genes involved in apoptosis, such as Bax, Bak, Bmf [24][25][26], and other genes involved in cellular signaling [27][28][29][30][31][32][33][34][35].
ADGRA3 (also known as GPR125; Fig. 1a) is an orphan member of the adhesion G protein-coupled receptor (ADGR/aGPCR) superfamily, belonging to the ADGRA subfamily (consisting of members ADGRA1-3) [36].ADGRA3 was initially described as a marker for spermatogonial stem cells [37][38][39], and we recently identified it as a factor of male infertility due to ejaculatory duct obstruction in half of all male mice lacking ADGRA3 [40].Moreover, Spina et al. recently described ADGRA3 expression on progenitor cells at the leading tips of the lacrimal ducts, which was required for proper tear film in mice, and on migrating progenitor cells at the invasive tips of ducts and branches during mammary development in mice [41,42].The direct cellular signaling mechanism for ADGRA3 has yet to be unraveled, but the receptor is expressed at the cell surface, where it constitutively internalizes [43].ADGRA3 interacts directly through PDZ interactions with both Disheveled (DVL) and Discs large (DLG) [44][45][46][47][48], indicating potential roles in both canonical WNT and non-canonical planar cell polarity (PCP) pathways.Supporting this, previous studies have shown that ADGRA3 negatively regulates the canonical WNT/β-catenin pathway [49,50], a pathway central to many developmental processes [51].
In this study, we aimed to investigate the role of ADGRA3 in female urogenital development and fertility.We hypothesized a similar role in females based on its role in the male reproductive system in mice [40] and its expression in spermatogonial stem cells in men [37][38][39].We, therefore, disrupted endogenous Adgra3 and found that 44% of ADGRA3-deficient females are infertile due to a lack of vaginal canalization at sexual maturity.Our data demonstrates a central role of ADGRA3 in distal vaginal tissue remodeling and formation of the vaginal opening.

Adgra3 is highly expressed in the prepubertal vagina and adult corpus luteum
To determine Adgra3 expression in the normal female urogenital system, we performed qPCR on tissues from 8-to 10-week-old C57BL/6J female mice in proestrus.Adgra3 was expressed ubiquitously in the urogenital system and endocrine tissues, including the pituitary and adrenal glands, with kidneys and ovaries expressing the highest levels and vaginae the lowest (Fig. 1b).In contrast, the other two ADGRA subclass members (Adgra1 and 2), though present in all tissues examined, were expressed principally in the pituitary gland and bladder (Adgra1) and to a lesser extent in the kidneys and ovaries (Adgra2; Additional file 1: Fig. S1a and b).Adgra2 expression was similar to Adgra3 in most other urogenital tissues (Fig. 1b and Additional file 1: Fig. S1b).Intriguingly, Adgra3 was expressed to a higher extent in the ovaries (1.6-fold) and vagina (7-fold) of prepubertal mice (postnatal day 15, P15) compared to 8-to 10-week-old C57BL/6J mice (p = 0.0001 and p < 0.0001, respectively; Fig. 1c).Given the high ADGRA3 expression in spermatogonial stem cells [37][38][39][40], we investigated ADGRA3 within the ovaries utilizing transgenic mice expressing ADGRA3 fused to β-galactosidase at the first transmembrane domain [38,[40][41][42].In these mice, ADGRA3 was predominantly located in the corpus luteum (Fig. 1d), with high expression in granulosa lutein cells (Fig. 1e and f ).
Fig. 1 Adgra3 is expressed during stages of vaginal remodeling and in the corpus luteum.a Illustrative structure of ADGRA3 (also known as GPR125).7TM, seven transmembrane domain; GAIN, GPCR auto-proteolysis inducing domain; HBD, hormone binding domain; Ig-like, immunoglobulin-like; LRR, leucine-rich repeats; PBM, PDZ-binding motif.b, c Endogenous Adgra3 expression in b whole tissues of the adult female urogenital system harvested from 8-to 10-week-old C57BL/6J females in proestrus (n = 4) and c in the ovary, pituitary gland, uterus, and vagina of premature (P15, n = 8) and mature (8-10 weeks, n = 8) Adgra3 +/+ mice in proestrus.Relative expression was calculated relative to the mean expression of housekeeping genes 36b4 and Ywhaz using the 2 −∆CT method, and significance was determined by unpaired t-test (ovary: ***p = 0.0001; pituitary gland: ***p = 0.0003; uterus: **p = 0.0024; vagina: ****p < 0.0001).Each dot represents an individual animal.Data are presented as mean ± SEM. d ADGRA3-β-gal expression in whole ovarian mount from an adult Adgra3 lz/+ mouse.The asterisk (*) indicates the ovary.U, uterine horn.e, f Representative images of a paraffin-embedded section of X-gal-stained ovarian whole mount from an Adgra3 lz/+ mouse (n = 3), with expression located in the granulosa lutein cells within the corpus luteum (f).CL, corpus luteum This phenotype was not observed in Adgra3 +/+ (0/58) or Adgra3 +/− female littermates (0/100) (Additional file 2: Fig. 2a).Obstruction of the vaginal entrance resulted in intrauterine fluid accumulation, distention of the uterine horns (hydrometrocolpos) and vagina (hydrocolpos) (Fig. 2e), and gross enlargement of the female reproductive tract (Fig. 2f ).Serial histological sectioning confirmed complete closure of the vaginal entrance in Adgra3 −/−CV females (Fig. 2g).The rest of the reproductive tract developed normally in all genotypes with no macroscopic or histological abnormalities.However, the pressure in the lumen distorted the vaginal and uterine histology of Adgra3 −/−CV females (Fig. 2g and Additional file 3: Fig. S3a-b).Moreover, the body weight of 8-to 10-week-old Adgra3 −/−CV females tended to be slightly higher, likely due to intrauterine fluid accumulation within the uterine horns, as the weight of other organs, such as the ovary, kidney, and brain, did not differ between genotypes (Additional file 3: Fig. S3c-g).

Ovarian morphology, estrus cycle, and puberty onset are unaltered despite ADGRA3-deficiency
In mice, vaginal opening is a hallmark of puberty, triggered by the first ovarian cycle and the corresponding increase in estrogen [11,12].As a measure of ovarian function, the onset of puberty was determined by the first day of the vaginal opening, which was similar between Adgra3 −/− (n = 5) mice and Adgra3 +/+ (n = 14) littermates (28.2 ± 0.66 and 27.6 ± 0.67 days, p = 0.6480).Furthermore, the length of the estrus cycle and the frequency of each phase of the estrus cycle did not differ between Adgra3 −/− and Adgra3 +/+ littermates (Additional file 3: Fig. S3h-i).These studies could not be performed on the infertile Adgra3 −/−CV females with closed vaginae.Nevertheless, by comparing the ovarian morphology of both Adgra3 −/− and Adgra3 −/−CV mice to that of Adgra3 +/+ mice, we found no difference, as follicles in all stages, including the corpus luteum, were observed in ovaries from all three groups (Fig. 3a).

Estrogen treatment fails to induce vaginal opening in mice lacking ADGRA3
To exclude reduced estrogen as the cause of the imperforate vagina, the impact of exogenous estrogen treatment was investigated in pre-pubertal mice.Vehicle-injected Adgra3 +/+ and Adgra3 −/− mice at P15 maintained a closed vagina (Fig. 3e), whereas all Adgra3 +/+ had a fully formed vaginal entrance after 3 consecutive days of estrogen administration (Fig. 3f ).However, the vagina remained closed in 55% (6/11) of Adgra3 −/− females (Fig. 3f ), suggesting that the phenotype was not due to low endogenous estradiol production at puberty onset (55% vs. 44% in the original colony, p = 0.4811, χ 2 ), but rather to a lack of estrogen responsiveness, even when circulating estrogen levels were high.In support of this, we observed significantly lower femoral bone mineral density (BMD), a long-term proxy of estrogen and estrogen signaling, in 1-year-old female Adgra3 −/− mice compared to Adgra3 +/+ mice (0.054 ± 0.0014 vs. 0.060 ± 0.0010 g/cm2, p = 0.0198; Fig. 3g).
As the development of the female reproductive tract is a complex process with many contributing factors, we measured an array of steroid hormones and compared the expression patterns of relevant mRNA in prepubertal (P15) and mature (8-to 10-week-old) female mice (Additional files 4, 5 and 6: Tables S1-S3).No differences were found between genotypes regarding androgen levels (Additional file 5: Table S2) or expression of androgen receptor (Ar) or anti-Müllerian hormone receptor type 2 (Amhr2) at P15 or 8 to 10 weeks  S1 and Additional file 6: Table S3).The anogenital distance (AGD), a measure of adequate androgen exposure, did not differ between Adgra3 −/− , Adgra3 −/−CV , and Adgra3 +/+ mice at any age (Additional file 7: Fig. S4).

BCL2 family genes are upregulated in vaginal tissue from mice with imperforate vagina
Progesterone promotes epithelial apoptosis in the vagina [53].As the intrinsic apoptosis pathway is controlled by the B-cell lymphoma 2 (BCL2) family of proteins (Fig. 4e) [54,55], we investigated the expression of selected BCL2 family genes.We found a marked upregulation of Bcl2 in the vaginal tissue of mice with closed vagina (approximately 5-fold higher in 8-to 10-week-old Adgra3 −/−CV ) compared to Adgra3 +/+ mice (p = 0.0009; Fig. 4f ).Concurrently, higher expression of genes encoding apoptotic sensitizers, Bid (5-fold) and Bmf (3-fold), was found in vaginal tissue from Adgra3 −/−CV mice compared with Adgra3 +/+ mice (p = 0.0003 and p = 0.0129, respectively; Fig. 4g,  h).Importantly, the expression of all BCL2 family genes measured in vaginal tissue was not different between Adgra3 −/− and Adgra3 +/+ mice (Fig. 4f-h and Additional file 4: Table S1).Furthermore, before pubertal e Illustration of the intrinsic apoptotic pathway regulated by the BCL2 family.Extracellular apoptotic stimuli activate sensitizer proteins, including BID, BIM, PUMA, and BMF.The sensitizers either inhibit the group of apoptotic protector proteins (e.g., BCL2) or activate apoptotic executors (e.g., BAX and BAK1).Activating apoptotic executor proteins leads to the release of cytochrome c from the mitochondria and subsequent apoptosis [54].f-h mRNA expression of BCL2 family genes fBcl2, gBid, and hBmf in vaginal tissue from 8-to 10-week-old Adgra3 −/− (n = 4) and Adgra3 −/−CV (n = 3) mice relative to Adgra3 +/+ (n = 7) mice.Relative expression was calculated using housekeeping genes 36b4 and Ywhaz and the 2 −∆∆CT method.***p = 0.0009 (f), ***p = 0.0003 (g), p = 0.0129 (h).Each dot represents an individual animal.Data are presented as mean ± SEM.Significance was calculated using one-way ANOVA post hoc Dunnett's multiple comparisons test onset (P15), Bcl2, Bid, Bmf, Bax, and Bak1 were all expressed similarly in the vaginal tissue of Adgra3 −/− and Adgra3 +/+ mice (Additional file 6: Table S3).Other timed postnatal remodeling processes, such as eyeopening and release of the ear flap, were normal in all mice.No interdigital webs were observed in any mice during this study.Thus, changes in apoptosis factors correlated temporally with the manifestation of the closed vagina phenotype.

Discussion
Congenital abnormalities of the female reproductive tract are common and have severe consequences, but their etiologies are often unknown [7].Here, we report that ADGRA3 significantly influences female fertility and is crucial for normal vaginal canalization.Approximately half of the ADGRA3-deficient females failed to develop a vaginal opening during puberty, a phenotype found rarely (0.7%, 1/154) in wild-type C57BL/6 females [56].Adgra3 is expressed throughout the female urogenital system, with higher expression in pre-pubertal mice, which drops precipitously after sexual maturation.This is consistent with patterns of Adgra3 expression in the lacrimal and mammary glands, which also undergo significant postnatal development, where expression is associated with undifferentiated cell states and migratory processes associated with concurrent lumen formation [41,42].
In mice, vaginal opening correlates with a pubertal increase in plasma estradiol [12].However, we found no significant difference between adult Adgra3 −/− mice and littermate controls regarding puberty onset, estrus cycling, plasma estradiol, ovarian morphology, or Cyp19a1 expression.Moreover, exogenous estradiol treatment failed to rescue the vaginal phenotype in Adgra3 −/−CV females.Collectively, these results exclude reduced estrogen levels at puberty as the cause of the imperforate vagina.Furthermore, vaginal expression of estrogen receptors (Esr1, Esr2, and Gpr30) remained unchanged, indicating that the vaginal tissue's ability to receive estrogen signals remained intact in Adgra3 −/−CV mice.Nevertheless, ADGRA3 deficiency was associated with a significant reduction in established hallmarks of long-term estrogen signaling, such as femoral BMD, consistent with previous studies [40].However, as ADGRA3 was recently described as a positive regulator of osteoclastogenesis [57], and normal osteoclast activity is needed for dynamic bone remodeling, it is possible that the impaired BMD was also affected directly by the low ADGRA3 expression.
Estrogen has been linked to progesterone signaling in the uterus, ovaries, and mammary glands [58,59] and is a known inducer of Pgr expression in the vaginal epithelium and stroma [52,53].The marked reduction in vaginal Pgr mRNA expression in adult Adgra3-knockout mice with closed vaginae (Adgra3 −/−CV ) but not in those with open vaginae or in wild-type littermate controls supports an overall impairment of the downstream effectors of the estrogen axis.Notably, the highest expression of Adgra3 was found in the ovary, specifically in the corpus luteum, the location of progesterone production [60].As the primary role of progesterone is to support the embryo and secure implantation in the uterus [60], a suboptimal progesterone axis in Adgra3 −/− dams during gestation could explain their smaller litter size.In support of a connection between ADGRA3 and progesterone signaling, the ocular phenotype observed in Adgra3 cre/ cre and Adgra3 lz/lz mice exacerbates in dams during pregnancy when progesterone levels are high but resolves during weaning [41].However, downregulation of sex hormone receptors is not sufficient to drive an impairment of vaginal canalization, as evident from no reports of imperforate vagina in single-and double-knockout of Esr1 and Esr2 [61] or Pgr-null mice, though these mice have severe fertility issues [58].This suggests that other factors, driven by the lack of ADGRA3, contribute to the observed vaginal phenotype in our mice.
Although the exact mechanism of estrogen-driven apoptosis during vaginal canalization is unknown, progesterone signaling has been linked to apoptosis in the vaginal epithelium [53].We observed no significant difference in apoptosis factors during the prepubertal period but found an overall increase in vaginal Bcl2, Bid, and Bmf expression in 8-to 10-week-old Adgra3 −/−CV mice.Due to the accumulated fluid within the reproductive tract, the upregulation of Bcl2, Bid, and Bmf may have been secondary to the pressure within the lumen rather than a direct loss of Adgra3.However, a 1.4-fold increase in vaginal Bcl2 expression was observed in Adgra3 −/− at P15 (Additional file 6: Table S3).Although non-significant, this timing supports the concept that ADGRA3 may function to regulate the balance between pro-and anti-apoptotic signals during vaginal canalization.Several apoptotic genes are associated with vaginal aberrations in mice.For example, triple-knockout of apoptotic sensitizers Bid, Bim, and Puma or dual knockout of pro-apoptotic genes Bak and Bax causes imperforate vagina in 100% of knockout mice, and this phenotype also occurs with lower penetrance in 20% of Bmf-knockout mice [24][25][26].Furthermore, 100% of transgenic mice overexpressing the anti-apoptotic gene Bcl2 have a closed vagina in adulthood and, similar to our findings in Adgra3 −/−CV mice, cannot be rescued by estrogen treatment [11].Thus, in contrast to targeted deletion of the sex hormone receptors, which is insufficient to drive imperforate vagina development, changes in apoptotic regulators are strong enough to drive this phenotype.
Development of the Müllerian ducts and female reproductive tract is highly dependent on proper WNT signaling [9,15], as evident from imperforate vagina having a reported penetrance of 30% to 100% in knockout experiments of other genes involved in WNT signaling (e.g., Ctnnb1, Wnt4 and Wnt5a, Dlg1, Vangl2, Map3k1, Ovo1) [16][17][18][19][20][21].Notably, a missense mutation in β-catenin results in a similar phenotype in male mice, as we previously described for Adgra3 −/− males [40], in addition to the imperforate vagina in female mice [16].At the signaling level, it is still unknown whether ADGRA3 signals via G proteins.However, several groups have reported a strong linkage between ADGRA3 and WNT signaling, focusing on both the developmental and other roles of this receptor [44][45][46][47][48][49][50].A similar linkage with WNT has also been reported for the closely related ADGRA2 [62].Similar to ADGRA3, ADGRA2 is required for tip cell function, although in endothelial cells during angiogenesis in the brain [62].Intriguingly, vaginal Adgra2 expression is slightly higher in P15 Adgra3 −/− mice and significantly higher in 8-to 10-week-old Adgra3 −/−CV mice compared to littermate controls (Additional file 4: Table S1 and Additional file 6: Table S3).Further investigation is required to determine if ADGRA2 compensates for the loss of ADGRA3 and if this is the reason for the incomplete penetrance of the vaginal phenotype.

Conclusions
Our findings suggest that ADGRA3-deficiency impairs the estradiol-dependent vaginal canalization via an altered balance in apoptotic regulators, potentially driven by WNT signaling.Future efforts will be directed at deciphering if and how ADGRA3 and ADGRA2 act in concert to affect WNT-driven epithelial remodeling in the distal vagina.As G protein-coupled receptors are excellent drug targets, our study provides support for investigating ADGRA3 as a target in diagnostic screening and/or therapy for obstructive vaginal malformations in humans.

Animals and housing
Adgra3 lz/+ mice were described by Seandel et al. [38] and Spina et al. [41,42].Adgra3 −/− and Adgra3 +/+ mice were generated from heterozygous breeding and were used previously in Nybo et al. [40].All mice were housed in individually ventilated cages with 8-10 air changes per hour.The environment was maintained at 22 °C (± 2 °C) and 55% (± 10%) air humidity.The stables maintained a 12:12-h light-to-dark cycle.The mice received an ad libitum chow diet (D30, SAFE, Rosenberg, Germany) and had free access to tap water.Housing, breeding, and experiments were conducted according to institutional guidelines and approved by the Animal Experiments Inspectorate under the Danish Ministry of Food, Agriculture, and Fisheries.License numbers for this study were 2017-15-0202-00117, 2017-15-0201-01235, and 2018-15-0201-01442.

Genotyping
Genotypes were determined using genomic DNA from ear or tail biopsies.Tissue samples were lysed in boiling (90-95 °C) 25 mM sodium hydroxide (NaOH) and 0.2 mM Ethylenediaminetetraacetic acid (EDTA) solution for 30 min, and DNA stabilized in 40 mM Tris-Hydrochloride.Three primers (forward P2: CTG AAG TTC TTT TGA CAG AAT CTC GGCAC, forward P5: AGA ACC ATG AAA CGT GGG ATA CCT GTTTC, reverse P3: TGG AGC ATA CAC GAG CAC TCT GTT A GTCA) were used to identify the wild-type (200 bp) and mutated alleles (500 bp).PCR was performed using PCR Master Mix (2X) (Thermo Fisher Scientific, Vilnius, Lithuania) according to the manufacturer's protocol.Primer products were amplified using an Eppendorf Mastercycler EP Gradient S Thermal Cycler with initial denaturing at 94 °C for 2 min, followed by 30 cycles of 94 °C for 30 s, 65 °C for 30 s, and 68 °C for 2 min, and completed at 68 °C for 8 min.PCR products were loaded into 1% agarose gel containing 1× GelRed Nucleic Acid stain (10000X) (Biotium, Fremont, CA, USA) alongside 1 kb Plus DNA Ladder (Thermo Fisher Scientific, Vilnius, Lithuania) and separated by electrophoresis.Products were detected in the gel using 1× Gel Loading Dye Purple (6X) (New England BioLabs, MA, USA), and PCR products were visualized using a Protein Simple SA-1000 (Red) gel imager.

Fertility trial
The fertility of female mice was studied by housing 8-week-old Adgra3 −/− and Adgra3 +/+ mice together with 8-week-old wild-type males for 5 months.Breeding pairs of four (one female Adgra3 +/+ , two female Adgra3 −/− , and one male Adgra3 +/+ ) were caged together.Adgra3 +/+ females were housed with males 5 days ahead of Adgra3 −/− females to prevent simultaneous birthing.Litter sizes were monitored at birth during the study period.

Monitoring puberty onset and estrus cycle
To determine the onset of puberty, the status of the vaginal opening was evaluated each morning from P20 until complete formation or P35.The estrus cycle stage was ascertained for 17 consecutive days from the age of 7 weeks by vaginal lavage and vaginal cytology as described previously [63].Briefly, the vagina was flushed with 100 μl saline, and the smear was air-dried, stained for 1 min in 0.02% crystal violet in absolute ethanol, and rinsed twice for 1 min in distilled water.The stages of the estrus cycle were defined by the presence of primarily anucleate cornified cells (estrus), a combination of anucleate cornified cells and leukocytes (metestrus), mostly leukocytes (diestrus), or the presence of primarily small, nucleated cells (proestrus).Vaginal lavage was performed in the morning between 9 and 11 am and in the evening between 8 and 10 pm, as the metestrus and proestrus stages of the cycle last less than 24 h in mice [63].

Pre-pubertal estradiol treatment
In mice, canalization of the vagina occurs at puberty via an estrogen-mediated induction of apoptosis of the distal vaginal epithelium to generate the vaginal entrance [10,11].Thus, the opening of the murine vagina remains closed at birth and is sensitive to local and systemic alterations in estrogen signaling that accompany sexual maturation [12][13][14].Pre-pubertal mice were injected intraperitoneally with 2.5 μg/g 17β-estradiol (Sigma-Aldrich.St. Louis, MO, USA) dissolved in isopropanol and diluted in flax seed oil for three consecutive mornings starting from postnatal day 12 (P12).The status of the vaginal opening was assessed and recorded each morning.

Anogenital distance (AGD)
The AGD of female Adgra3 −/− mice and control littermates was assessed by photograph, allowing for repeated assessment.The mice were placed on a Plexiglas platform with a built-in measuring tape and a digital camera approximately 10 cm below.The mice were held stationary by the tail and photographed.These photographs were depicted at high magnification and analyzed by two blinded investigators.The AGD was measured as the length of the perineum from the center of the genital papilla to the proximal end of the rectum [64].

Blood sampling, tissue preparation, and qPCR
Adult animals were euthanized by cervical dislocation during the proestrus stage.Blood was drawn from the main aorta immediately after decapitation, and organs were harvested and snap-frozen in liquid nitrogen or immersed in 3.7% formaldehyde in phosphate-buffered saline (PBS).Plasma was separated from whole blood by centrifugation at 1000 rpm at 4 °C for 15 min.The right femur was cut free from the muscles, fixated in 70% ice-cold ethanol, and stored at 4 °C until further analysis.Total RNA was extracted from snap-frozen tissue using the RNeasy Mini Kit (Qiagen, Hilden, Germany) or RNeasy Micro Kit (Qiagen, Hilden, Germany) for tissues > 5 mg or < 5 mg, respectively, according to the manufacturer's protocol.RNA was measured using NanoDrop One (Thermo Fisher Scientific, Wilmington, DE, USA), and reverse transcription of purified RNA (2000 ng) was performed using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA) according to the manufacturer's instructions.Real-time quantitative PCR (qPCR) analyses were performed using the Quantstudio 6 Flex Real-Time PCR System (Applied Biosystems, Forster City, CA, USA) and Precision-PLUS qPCR Master Mix with low ROX and SYBRGreen (Primer Design, Chandler's Ford, UK).All primers were used at a concentration of 250 nM per 10 μl reaction, and cDNA was used at 1 ng/μl.Gene expression was reported relative to housekeeping genes tyrosine 3/tryptophan 5-monooxygenase activation protein, zeta polypeptide (Ywhaz), and ribosomal protein, large, P0 (36b4/Rplp0), and was calculated by the Livak (2 −ΔΔCt ) method (for two or more groups, normalized to Adgra3 +/+ ) or the 2 −ΔCt method (for one group).Primers were ordered from TAG Copenhagen (Frederiksberg, Copenhagen), and their efficiencies were established according to MIQE guidelines [65].Primers with efficiency (E%) between 90 and 110% were used.A complete list of primer details is provided in Additional file 8: Table S4.

Dual-energy X-ray absorptiometry (DEXA)
Femora from 1-year-old Adgra3 +/+ and Adgra3 −/− females underwent densitometry using a DEXA scanner (Sabre XL, Norland Stratec, Pforzheim, Germany).The entire femur was scanned at a pixel size of 0.1 × 0.1 mm and a velocity of 3.0 m/s.The BMD was determined using the software supplied with the scanner.Quality control was conducted using two solid-state phantoms.

Hematoxylin and eosin staining
Formalin-fixed tissues were paraffin-embedded, sectioned (4-to 5-μm), mounted on glass slides, and stained at the Histolab, Core Facility of Integrated Microscopy (Biomedical Institute, University of Copenhagen, Denmark).Hematoxylin and eosin staining was used to determine organ and tissue morphology.Serial sections were used to determine the nature and extent of vaginal obstruction.Images were acquired using an Axio Scan.Z1 (brightfield) Slide Scanner (Carl Zeiss, Jena, Germany) and Zeiss Zen Blue 3.0 Software (Carl Zeiss, Jena, Germany).

Statistical analysis
All statistical analyses were performed in Prism 9 (Graph Pad Inc., San Diego, CA, USA).The qPCR data (values from cycle threshold) were normalized to internal controls Ywhaz and 36b4.Outliers in qPCR analyses were identified using the ROUT method with a false discovery rate of Q = 1%.Outliers were excluded from further analyses and figures.Calculations were based on ΔCT values for qPCR data.Before statistical analysis, the normality of data distribution and homogeneity of variances were assessed using Shapiro-Wilk tests and Levene or Bartlett tests, respectively.Variance and p-values were calculated using a two-tailed, unpaired Student t-test (two groups) or one-way ANOVA followed by Dunnett's multiple comparisons test for post hoc analysis (three groups).When the sample size was too small to assess normal distribution, the nonparametric test of Kruskal-Wallis, followed by Dunn's multiple comparisons test, was performed to provide a mean of statistical analysis without assuming a specific distribution.Observed vs. expected characteristics (observed closed vagina with vs. without estradiol treatment) were compared using the χ 2 test.Continuous outcomes are presented as mean ± SEM, and significance was defined as p ≤ 0.05.

Fig. 2
Fig. 2 Infertility in Adgra3 −/− females is due to an imperforate vagina.a Percentage of fertile Adgra3 +/+ (n = 4) and Adgra3 −/− (n = 9) females.b Number of litters per Adgra3 +/+ and fertile Adgra3 −/− mothers after 5 months of consecutive breeding.Each dot represents an individual animal.c Number of pups in each litter from fertile animals (*p = 0.0312).Each dot represents one litter.Data are presented as mean ± SEM, and significance was determined by an unpaired t-test.d, e Representative images of vaginal openings (d) and female reproductive tracts (e) in 10-week-old Adgra3 +/+ , Adgra3 −/− , and Adgra3 −/−CV mice (n = 3 per group).Arrow indicates an opened vaginal opening, and an arrowhead indicates a lack of vaginal opening observed in 44% of all Adgra3 −/− females in the colony.An asterisk (*) indicates the ovary.U, uterine horns; V, vagina.f Dissected female reproductive tract (from one ovary to the vagina) of 10-week-old Adgra3 +/+ and Adgra3 −/− mice compared to the enlarged reproductive tract of a 10-week-old Adgra3 −/−CV .g Hematoxylin and eosin (HE)-stained sections of the distal vagina in 10-week-old mice.Animals with an open vagina were euthanized after a smear test confirmed the proestrus stage